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FOREWORD 

This report was prepared in fulfillment of the requirements of NASA 
contract NAS8-29193- George C. Marshall Space Flight Center was the con- 
tracting agency. Technical Coordination was maintained through Dr. Steve 
Winder, Chief, Statistical Dynamics Section, S&E-AERO-DDD. 



NORTHROP SERVICES, INC. 


TR-1246 


ABSTRACT 

An investigation was made of the problems of excess control devices and 
insufficient trim control capability on shuttle ascent vehicles. The trim 
problem is solved at all time points of interest using Lagrangian multipliers 
and a Simplex based iterative algorithm developed as a result of the study. 
This algorithm has the capability to solve any bounded linear problem with 
physically realizable constraints and minimize any piecewise differentiable 
cost function. Both solution methods also automatically distribute the com- 
mand torques to the control devices. 

It is shown that trim requirements are unrealizable if only the orbiter 
engines and the aerodynamic surfaces are used. On the other hand if the solid 
engines are controllable there is ample control margin and realistic cost 
functions can be minimized to optimize the solution. 
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Section I 
INTRODUCTION 

This report presents the results of a study to develop a technique to 
solve the trim problem of an aerospace vehicle when the number of control 
variables exceeds the number of variables to control trim. Distribution of 
the trim commands to the various control variables in some "optimum" manner 
is a further goal of the study. The goal of the study is not a computer 
program as much as it is the development a practical method of understanding 
and solving the trim problem. 

The importance of solving this type problem became evident as the Space 
Shuttle vehicle requirements were established. In particular, the ascent phase 
of flight is one area in which these problems arise. In ascent the mated 
configuration has only one plane of symmetry causing the inherent problem of 
zero lift and zero moment occurring at different angles -of-at tack. Cross- 
coupling problems are also severe in some candidate designs. 

Several candidate configurations have insufficient control authority to 
control the vehicle in the presence of headwinds or crosswinds with orbiter 
engines only. When the aerodynamic surfaces on the orbiter are used, the 
shuttle may or may not be controllable, but in any case the number of control 
effectors becomes greater than the number of states to control, creating an 
infinity of trim solutions for any constant disturbance. If controllability 
problems require vectoring the solid-rocket motors, the number of trim 
combinations proliferate. 

The specific goal of this study was to solve the problem of distributing 
the trim commands to the various controllers, simultaneously using the extra 
control choice latitude to optimize some phase of the ascent flight. Further 
goals were that the solution to the trim problem be simple and easily imple- 
mented and that the solution process give the engineer a "feel" for the rela- 
tive difficulty of meeting the trim requirements as compared to the trim power 
available . 
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Of equal importance to solving this problem is the delineation of the 
concepts involved in its solution. Since the sample problems chosen are 
realistic, a computer program was necessary to handle the data efficiently, 
but the ideas behind the program receive the emphasis throughout this report. 


The problems that arise in trimming the Shuttle ascent vehicle are pri- 
marily in the roll/yaw planes. The shuttle lateral plane trim problem can be 
outlined in the following manner. 

• There are three vehicle state variables; roll, yaw, and sideslip 
angle. 

• There are eight control devices with effectiveness in the lateral 

plane: three orbiter engines, two solid-booster engines, right and 

left aileron and a rudder. Each engine has two degrees of freedom 
(pitch and yaw); however, the top orbiter engine is on the vertical 
centerline of the vehicle (049 data). Therefore, it has no effect on 
the lateral dynamics when it is pitched. The net result is 12 control 
states in the lateral plane. 

• All control devices are limited by control stops or hinge moments or 
both. 

• The problem is to find the range of control effector settings that 
satisfy the trim requirements on the states in some optimum manner. 

The reason for selecting the trim problem to optimize the control devices is 
that Shuttle command torques and forces were currently unknown. Therefore, 
trimming an aerodynamic disturbance, such as a mean wind, is a convenient sub- 
stitute for a torque command profile. The results of this study will give not 
only a trim solution, optimum in some manner, but a mechanization of the dis- 
tribution of commands to the control devices to achieve a given flight state. 


Figure 1-1 illustrates the approach described in the previous paragraph. 
The T^ quantities represent the torque command profile. The optimization 
scheme should be independent of the torque commands so that the guidance design 
functions and control design functions can proceed somewhat independently of 
each other until system integration time. 

The general problem can be stated formally as follows: 
minimize 

w(6^) 
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Figure 1-1 . FEEDBACK DIAGRAM SHOWING OPTIMUM SCHEME 


subject to the trim constraints 

<(.,(6.) - T = 0 

1 J <=1 

i = 1, 2, ...» tn 
j = 1, 2, . . . , n 


and inequality constraints 

6 . < 6 . < 6 . 

3 . — 1 — j 
'^min ^ '^max 

where the represent a torque command of a disturbance vector. 


In the Shuttle trim problem considered, the t>e represented by 

a set of i linear equations and the are just constants. This means that 
the equation involving can be solved for i *^j'^ chosen 

beforehand. The boundaries of the problem can be explored by setting j-i con- 
trol settings on the bounds and calculating the remaining i control variables. 


Data for a Shuttle type vehicle were furnished by the COR. These data 

are presented in Appendix A. The configuration is similar to the 049 config- 

2 

uration with a 400-ft forward fin to reduce the lateral stability. No drag 
data for the aerodynamic surfaces were included, but the rudder drag coeffi- 
cient was estimated from wind tunnel data on later configurations. The aileron 
drag coefficient was assumed to be equal to the rudder for the purposes of this 
study. No aileron side force coefficient was included in the data either, 
but this turned out to be unnecessary for the problems solved. 
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The position limits for the rudder and aileron are shown in the data as 
well as the hinge moment limits for these surfaces at various, flight times. 

The engine gimbal stop limits were assumed to be 10 degrees circular displace~ 
ment for the orbxter engines and 10 degrees diagonally up or down for the 
solid- rocket motors. No data on engine hinge moment limits were available so 
they were not considered in the problem solution. However, the method of 
solution used would handle this type of restriction with a simple data change. 
The mean cross-wind was used as a disturbance for much of the study (Appendix 
A). 


Figure 1-2 shows the mean crosswind included roll /yaw torques (normalized 
by dividing by q sb) at times 25-90 seconds. As can be seen, the yaw torque 
peaks at 60-65 seconds and the roll torque peaks at about 70 seconds. However, 
dynamic pressure peaks at 80 seconds so the maximum controllability problem 
occurs around 75-80 seconds. It is impossible to predict the exact time in 
advance, since the control problem also depends on the individual control 
effectiveness and their proportional roll/yaw effectiveness. These parameters 
are functions of thrust, dynamic pressure and Mach number. 

The approach and results section presents in somewhat chronological 
order the approaches taken and the results achieved in the study. Specifics 
involving data, some of the mathematical developments, and the computer program 
developed are included in the Appendices. 
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Section II 

APPROACH AND RESULTS 

The trim problem was approached by first formulating the constraints and 
limits for the control effectors in a form suitable for analysis. Then, to 
give the investigator a feel for the problem, several optimized trim solutions 
were calculated for special cases. After sufficient feel was obtained, the 
optimization problem was solved using a quadratic cost function. This problem 
has been solved many times and it can be shown that any cost functional that 
can be represented by a quadratic cost function can be solved by Lagrangian 
multipliers, the technique used in this part of the study. Several sample 
problems were reduced by . this method with the objective of giving additional 
insight into the problem. Since the Lagrangian multiplier technique cannot 
directly handle inequality constraints, this approach does not guarantee a 
realizable solution even if one is available. Therefore, the next step in the 
approach was to develop a means of solving the trim problem with inequality 
constraints and with a general form cost functional. 

Several methods were available to solve bounded problems with nonquadratic 
cost functions. All of them, except in special cases, require iteration to 
reach their solution. It was decided that a solution algorithm based on the 
simplex algorithm would be conceptually the simplest to implement. The advan- 
tage of simplex in the solution of the trim problem is that it can be treated 
as a black box which simultaneously solves the optimization problem with more 
variables than equations and keeps the solution within the variable bounds. 

Several sample problems were solved by this method and the results eval- 
uated. Finally, a dynamic survey was performed to give a feel for the inter- 
action of the Shuttle dynamics with the control trim minimization scheme, 

2.1 EVALUATION OF TRIM BOUNDARIES 

The solution to the trim problem (zero torques and forces) is an indefinite 
set of control trim placements. In fact there are an infinity of solutions, 
since there are fewer equations than control variables. However, the unknown 
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control variables are limited in their displacements so the set of possible 
solutions can be regarded as a bounded, infinite set. 

Figure 2-1 shows some of the boundaries of the trim solution as calculated 
for 80 seconds flight time. This figure gives an indication of the relative 
trim torques available at 80 seconds flight time. As can be seen, each control 
device was placed individually on its limit and the other control trim settings 
were computed. Since the possible trim solutions of 11 control devices form a 
volume in 11-space it is impossible, to show graphically; but the points shown 
in Figure 2-1 represent points on the boundary of the 11 volume. Since engines 
two and three do not pitch in the solutions and since they yaw about 75 percent 
of engine one they are not shown on their limits. 

Figure 2-2 shows the minimized control trim settings when each control 
device was shut out individually. This shows that the vehicle has sufficient 
control power for trim with any one controller out as long as the SRMs have 
TVC. It also shows that the aileron has little effect on the trim requirements 
since deletion of the aileron does not appreciably change. 

2.2 TRIM PROBLEMS WITH QUADRATIC COST FUNCTIONS 

A quadratic cost functional is a function made up of the sum of the squares 
of the variables in the function. The partial derivatives of a quadratic cost 
function are linear making the minimization of the cost a linear problem. 

Several techniques can be used to solve such problems exactly but the most 
convenient one for this problem was Lagrangian multipliers. 

The Lagrangian multiplier technique will minimize quadratic cost functions, 
giving an exact solution to the trim problem optimized with respect to the 
cost. Penalty functions can be devised to minimize such things as drag due to 
zero surface deflection, thrust gimbal angle, and off-nominal translational 
accelerations. 

The penalty functions can be accounted for in the cost functional. Defi- 
nition of the cost functional depends upon the particular problem being solved. 
In gene:ral, this functional is a summation of various cost terms which is to 
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Figure 2-2. CONTROLLABILITY WITH CONTROL DEVICES ZEROED INDIVIDUALLY 
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be minimized subject to the constraints. The individual terms differ from each 

other qualitatively and must be weighted in accordance with the requirements 

of the investigation. An example cost functional, oj, is defined below: 

2 2 

0) = K^x(DRAG DUE TO ELEVON DEFLECTION) + K2X (THRUST GIMBAL ANGLE) 

+ K^x(RUDDER HINGE MOMENT)^ + K^x(VARIATION OF CONTROLLER FROM 
MEDIAN VALUE) ^ + K^x (OFF-NOMINAL TRANSLATIONAL ACCELERATIONS)^, 

where the K scale factors are selected so that the various terms have the 
i 

proper relative costs. The and K2 terms are included so that the "optimum" 
trim solution will not produce excess AV losses. The term prevents a trim 
solution which requires an excessively large hydraulic system for rudder 
deflections. The term prevents a trim solution where a control variable 
value is close to its maximum limit. The term prevents a solution where 
the resultant translational accelerations are greatly different from the nomi- 
nal values. However, the Investigator must bear in mind that regardless of 
the quadratic cost functions devised, it will always reduce to the form 

m « 

u) = I (a. 6 .) 

iti J J 

because the K terms can always be summed with each other if the i control 
variable appears in the cost functional more than once. In effect, any study 
of quadratic cost functionals, no matter how exotic they may be, will always 
amount to a variation of the a^ coefficients, which are just constants. 

Since it was desired to develop a generalized technique for solving trim 
minimization problems and not to perform a sensitivity study on quadratic cost 
functions, it was decided to set up and solve a few sample problems to gain 
further insight into the general trim problem. No attempt was made to repre- 
sent all possible penalty functions in the cost functional. 


The problem was reformulated and Lagrangian multipliers were used to solve 
in the following manner: 
minimize 


m 


0 .( 6 .) 


= I 

j=l 
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subject to 


■^1 = ^ 'j ' b = ° 
j=i J ^ 

i=lj2, .•»>n • 


m 


Let F = 0 ) + I 
i=l 


To minimize F, set 
9F 


36 . 
J 


= 0 


2 - ly m 


This will yield m equations and m + n unknowns since the are not 
known. Switching to matrix notations: 


9F 

361 



< * 

3F 

36 

n 

<N 

II 

_0 '• 4- 



r 


i 

. 

6 

n 



■<=11 

^12 ^Im 


^1 

?21 

^22’. 

' 

: 

* 

c . 

’C 


X 

L 

niTL, 


V J 


define such, that 

2 
I ® 

~1 




0 




m 


let 


C = 


<=11 

S2---<=lm 



'y 

Si 


{6} 

i 

: 

^'nJ 

i 

c 

niQ 



n 


explicitly 


( 1 ) 


2-6 



NORTHROP SERVICES, INC. 


TR-1246 


r 

h 

T 


A = 


and (f) = 

/ 

I — 1 

^ . . 

1 .2 
T 


< 

.2 

A 

1 


iU 


returning to Equation (1) we can solve for 6 



also 

(b = C6-T = 0 
^ c 

eliminating X 

{6} = [B] [C]'’^ {[C] [£2] [C]’^} {T^} 


( 2 ) 

(3) 

(4) 


We have solved for 6 in terms of the C matrix, the penalty matrix U, and the 

t 

torque commands . 

An immediate result of using three constraint equations (yaw moment = roll 
moment = side force = 0) is that no practical minimum solution was obtainable. 
That is, the control settings are greater than the bounds. The zero side force 
requirement was the cause of this. When the requirement was dropped the prob- 
lem reduced to the following: 


minimize 




+(^6«e3p>'+("7^e3y>' 


(5) 


subject to 


', = c. <5 + c ■SR + c a + c a + a 

*■ ^6a ® ''aR 6ely ely ae2p e2p 6e2y e2y 

+ G„ a + c a , - (-C, B) = 0 

e3p *-ae3y 


( 6 ) 
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,=C a+C 6R+C 6ely + C 

2 Sa "«ely , 


6e2p + C 


(Se2y 


6e2p 


iSe2y 


+ C 


6e3p 


6e3p + C 6e y - (-C 6) - 0 

^6e3y 6 


( 7 ) 


where 3 comes from the mean wind. The a^ coefficients can be set to any value 
so that varying them will represent objective functions that can minimxze con- 
trol deflections, hinge moments, thrust losses due to drag and engine displace- 
ments. This technique cannot automatically handle the limits on the control 
devices so an iteration procedure was devised to solve for 6 and then check for 
excessive values. If they occurred the penalty coefficients were systematically 
varied to see if the solution could be driven within the physical bounds. 


Using the 049 data, the overriding problem near max Q (75-80) seconds is 
controllability. By choosing the penalty coefficients according to the fol- 
lowing formula a solution that minimizes control deflections can be obtained. 

a. = 10/minimum (|5. |or|6 |) . 

^ ^max min 


The philosophy here is that the control devices that are most severely limited 
will be most heavily weighted. Unfortunately, the solution to this problem, 
in the neighborhood of max Q, yielded a rudder setting approximately two times 
as large as the control power can furnish (hinge moment limited) , The aileron 
was also at or near the limit at this time. Systematically varying the penalty 
coefficients to drive these two devices within the control limits resulted in 
exceeding the control limits on all of the controllers . While this is not 
formal proof that there is no solution at these points it is very strong evi- 
dence for this conclusion. 

Additional insight into the problem can be gained by showing that (for 
the 049 data) engines 2 and 3 are not independent of each other if there is 
no reason to weigh one's displacement more heavily than the other. 


From the data the following holds: 
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M 


Z’ T 


T 


r/ 


T_ + T 


11 c 


12 c 


J 


{hi\ * ^1%) 


a 

1^ 


and if and and (equal weights on engines 2 and 3). 


Equation (4) can be written as 


*ely 


e?.p 


e2y 


e3p 


’e3y 



r^c 1 



b)W H' 


i = 

fc. a + c y\ 

V. ^6a , ^6a ) 


^22 vS 


/c„ a + c y\ 
\ ^6R “<SR ) 


33(''£ 


ot + C 


6ely 


5ely 


') 


^ Se 2 p ^) 


^44( 

”55 


a - C 


6e2p 


n 


6e2p 


a + C 


6e2y 


n. 


6e2y 


) 

') 


( 8 ) 


6e2p = “6e3p = ^ 


6e2y = -6e3y = Q 


55 


( 


a + C Y 
6e2p ^6e2p 


a + C 


n 


6e2y 


By a similar approach it can be shown that engine 1 is not independent of 
engines 2 and 3 if its penalty coefficient is the same as 2 and 3. The 
following ratios hold: 

^ilely " ^£e2y (Ze2‘^Zcg)cos Pe2CosYe2 


(9) 
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cosP 


ei 


n 


6ely 


n 


6e2y 


cosPg2^°®^e2 


Uo) 


where Is the Z location of engine one, is the pitch angle of engine 

one and is the yaw angle of engine 2, The result is; 




cosP 


el 


cosP 


ei. 


ely 


(Zp«+Z„^) cosPeoCOsY 


' 62 -cg 


®2 *2 


a C + Y C 

0 cosPo^cosY - rix o 

5 e 2 y ^^=*^62 e2 6e2y 


6e2y 


'S.e2y 


+ C 


n 


6e2y 


J(ll) 


The relationship between 6ely and 6e2y is not a direct ratio because A and Y 
are functions of T^ and T^, which are the disturbance torques. If the relative 
value of the disturbance torques change with time, or a different kind of dis 
turbance is used, the ratio between 6ely and 6e2y will also change. 


Similarly it can be shown that the left and right ailerons operate as 
mirror images of each other. Since it is difficult to imagine a case where it 
would be advantageous to operate them independently they have not been sepa- 
rated in any of the work done. However, the orbiter engines have all been 
operated both in pairs and independently in optimization procedure to date 
and the results have verified the developments in the previous paragraphs. 


Figure 2-3 shows the results of minimizing the control trim displacements 
at selected time intervals between the flight times of 25 and 90 seconds. As 
can be seen, the trim settings exceeded the hinge moment limits on the rudder 
and, to a lesser degree, on the aileron near maximum dynamic pressure. The 
severity of the controllability problem is illustrated in Figure 2-A. Since 
the rudder and the aileron exceed the hinge moment limits at 75 seconds, this 
time point was chosen as a test of the penalty coefficients iteration scheme. 
If a control device, say 6^, exceeded its maximum absolute value then its new 
penalty coefficient (aj) was reset by the following formula; 



( 12 ) 


where min [ ] means the minimum of the quantity in brackets. 
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Figure 2-4 is a plot of the iteration steps the computer program followed 
at 75 seconds flight time. Iteration forced the rudder down near its maximum 
but also forced the engines to exceed their limits. The program failed to 
find a realizable solution, indicating the magnitude of the controllability 
problem at this time point (with only a mean wind disturbing function) . As 
previously stated, this does not prove that the trim tortjues are nonrealizable 
but points strongly to the possibility. 

To alleviate the controllability problem the solid rocket motors (SRM) 
were allowed to gimbal (TVC) . Again the engines acted as a pair as long as 
the penalty coefficients on the engines were equal. As Figure 2-5 shows, the 
solid engines relieved the controllability problem. The peak rudder deflec- 
tion of 3.0 degrees was only about 13 percent of the maximum rudder deflection 
(hinge moment limited) at that time. The peak solid engine deflection was 
about 27 percent of the maximum (10 degrees) at 70 and 80 seconds. 

The conclusions drawn from the Lagrangian approach to trim optimization 
are as follows : 

• Controllability is a problem on the 049 even with the aerodynamic 
surfaces, if SRM TVC is not used. 

• If SRB TVC is used controllability is no longer a problem. 

• The aileron is relatively ineffective in trimming the vehicle. 

• If there are no reasons to weight them differently, engines two and 
three yaw equally and pitch oppositely. The solid engines pair also 
with the right engine displaced right and down and left engine right 
and up . 

• Trimming the side force is probably impractical, 

I 

In summary, the Lagrangian approach gives insight into the problem and an 
exact solution to the quadratic cost problem. However, if the cost can not be 
adequately represented by a quadratic the solution will be useless. Further- 
more, the procedure does not guarantee a solution within the variable bounds 
nor does It Indicate if a feasible solution is possible. 


PRECEDING PAGE BLANK NOT FIT.MED 
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2.3 BOUNDED TRIM PROBLEMS WITH GENERAL COST FUNCTIONS 

The trim problems investigated in this study had control variable bounds 
that typically limited their solutions. In the first part of this study it 

was necessary to obtain a solution and then check to see if the solution was 
physically realizable, that is within the variable bounds. Furthermore if the 
actual cost function was not quadratic then the Lagrangian multiplier technique 
could not be used at all. 

A solution process analogous to the Lagrangian multiplier technique can be 
used to solve this class of problems . However the details of the solution are 
involved and the solution process does not provide insight into the problem, 

A brief description of this approach is included here for clarity. 

If the cost functional is differentiable, but the derivatives are non- 
linear, setting the partials of the cost with respect to the control variables 
to zero yields a set of nonlinear homogeneous equations which can be solved by 
iteration (if the partials can be written analytically) . If the partials can 
only be represented by graphs or tables the mechanics of the solution are much 
stickier . 

Handling the inequality constraints adds more complexity to the problem. 
This requires a transformation of variable and definition of n additional 
variables linearly related to the original n variables. 

Since one of the goals of this study was to develop solution concepts, it 
was felt this brute force approach would not totally fulfill the work statement. 
Therefore, the following linear programming technique utilizing simplex was 
devised. 

The idea behind this approach was that the simplex algorithm can be used 
to force an iteration procedure to converge in the vicinity of the lowest cost. 
During the convergence procedure simplex will also require the controller 
bounds to be observed. The cost functional in a simplex problem may be recog- 
nized as a summation of partial derivatives. Conceptually it is very simple 
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to place any form of partials in the simplex cost row and the simplex solution 
that results will be as good as the values of the partials over which the 
iteration ranges . If this range is held small enough the solution will approxi 
mate the minimum cost on that range. 


The primary difference between a linear programming approach to the prob- 
lem and the Lagrangian multipler technique is that the control limits are 
considered in the solution by linear programming. In other words linear pro- 
gramming will produce a feasible solution in proximity to the optimal solution 
If the control limits are high enough to not impact the solution both methods 
will give approximately the same answer. 

As stated previously the ultimate goal of this study was the development 
of concepts and techniques to solve the trim problem. The solution to the 
general trim problem is actually an algorithm developed around the simplex 
technique. Therefore this algorithm is the primary result of the study. The 
central philosophy of this algorithm is to use the simplex procedure (ref. 1) 
to solve the problem of minimizing a general cost function^ £ii, over a small 
range of the variables' permitted values. 

2.3.1 Description of Solution Algorithm 

The mathematical basis for the solution algorithm rests on linearizing 
the cost function, w(6^), about a feasible solution, The cost can be 

written as follows! 


i(6 ,) = o)(6^ ) + I 

^ o i=l 


dm 


( 6 . - 6 . ) 

i o 


The problem is to find a way to minimize m(6^) (which is always positive of 
course) . This problem is the same as making 


I 


i=l 


dCJ 

86 



6 . ) 
1 . 
o 


the largest negative value possible. As will be shown later. Simplex will 
solve the latter problem, but an iteration algorithm must be used to find the 
minimum oj(6^) because it is nonlinear in the general case. 
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The algorithm developed to solve the trim problem can be described in 
simple terms. The simplex algorithm itself is described in Appendix B, After 
the problem is defined the variables are transformed so that they are all non- 
negative, which is a requirement of simplex. Then, if the neighborhood of the 
solution is known beforehand, an approximate solution is used as a starting 
point. If nothing is known about the solution the entire range of the vari- 
ables can be allowed in the starting solution. In either case a feasible 
(meets the trim requirements and variable bounds) solution is calculated to 
start the procedure. Since simplex utilizes a cost row made up of linear costs, 
this row must be filled with the partial derivatives of the cost function with 
respect to the appropriate variables. The simplex procedure will solve for the 
minimum cost solution based on the linearized costs and keep the solution within 
the bounds of the variable limits. Of course the solution is only as good as 
the linear approximations to the cost partials. To converge on the correct 
solution temporary bounds, which are always within the real bounds, are created 
about the feasible solution. These bounds are set using a range percentage, 
delta, about the feasible solution. After the simplex solution is calculated, 
using the temporary bounds, the actual total cost, u), is computed. As long as 
the new total cost improves the range percentage remains the same. The problem 
is set up and solved again, using the last computed solution as a starting 
point about which the new bounds are set and at which the new linearized partial 
costs are computed. If the total cost does not improve delta is halved and the 
new solution is calculated. When the range percentage becomes small enough the 
simplex solution will be approximately the exact solution. At some improvement 
criteria the process is terminated . A simplified flow chart of the procedure 
is shox<m in Figure 2-6 . 


2.3.2 Formolotion of Simplex Problem 

The classical linear programming problem is as follows: 

Find > 0 such that 


Z = y a? is minimum and 

j=i ' 


(13) 

(14) 
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• i -t 

I ^ b (i = 1, m; m < n) 

j=l J 


(15) 


The are our control trim settings, Z is the objective function and Equation 
(15) represents the zero torque requirements and the control limits. 


The formal statement of the linear programming problem indicates that 
only positive values of the problem variables are allowed. Furthermore, the 
cost function, Z, is a linear summation of the products of the problem vari- 
ables and their constant multipliers. Since any practical control trim problem 
may require some negative control settings, a transformation of variables was 
performed on the problem variables. 


If we choose X. = 6. - 6, 
3 3 


^mln 


then since 


6 . < 6 . < 6 ^ 

^min ^ ^ 


max 
- 6 , 


0 ^ X. < 6 

^ ^max ^min 


and since 


6 . ^ 6 . ^0 
^max ; ^min 

the requirement on X^ is consistent. 


(16) 


In the notation of Equations (13) through (15) the problem has m con- 
straints and n variables. From experience the algorithm will solve the problem 
in ra to 3m steps (ref. 2). However, in the trim optimization problem simplex 
has to be applied iteratively until the cost reaches a near minimum. If the 
original solution was allowed to range over the whole boundary, the investi- 
gator's experience showed simplex would be applied as many as 20 times. If 
the solution is approximately known beforehand, the number of iterations was 
reduced sharply. One way to predict the approximate solution is to use the 
solution from the previous nearby time steps. 


The cost functional can be any piecewise differentiable function. A 
description of the computer program used to solve the trim problem is in 
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Appendix C. The cost functional is represented by two simple subroutines so 
that if it needs to be changed in form (not just coefficients) it can be 
readily reprogrammed. 

For purposes of testing the validity of this approach as well as obtaining 
useful information, a cost function was defined to represent the drag losses 
due to aerodynamic surface deflections and the thrust losses due to vectoring 
the engines off the centerline. The drag losses are represented as typical 
linear aerodynamic terms which are functions of the absolute aero surface 
displacement. The thrust losses are calculated by subtracting the cosine of 
each engines' circular displacement from unity and multiplying by its thrust. 
Since each engine is represented in the control equations by a pitch and a yaw, 
the circular displacement is the square root of the sum of the squares of the 
pitch and yaw displacements. The linear programming procedure will handle 
any piecewise differentiable cost function but this cost function in Equation 
(16) was chosen in the bulk of the study 


i=l 6 i=3 , 5 , 7 , 9 » , »n 


The simplex costs are defined as 


Z = 




36 . 

X 


6 . 

X 


solution 


or 


Z = [qS^^. C 


ref 1\T ' i 




sim;o_. + 6. 

2 2 


= 3, 5, 7, 9 ... 

(17) 


sin 


+ T. 


+ -S? 1 

» X x-1 


2 2 
'i ^ 'i-1 


6. i = 4,6,8,10, 
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The rudder drag coefficient was set to 0.7 x 10 nt/deg, a value that was 
estimated from looking at Shuttle rudder data for a similar configuration. 

For convenience, since there was no aileron drag data available, the aileron 
drag coefficient was set to that value also. 

2.3.3 Minimum Thrust/Drog Losses Withoug SRM TVC 

Figure 2-7 shows the minimum thrust, minimum drag results for 50 to 90 
seconds for the 049 configuration w The disturbance is the same as that used 

V ..k ^ 

before, the mean wind. This plot shows the trim results when the solid engines 
are not used. The "infeasible" region is the time between 74 and 85 seconds. 
During that period trim requirements cannot be fulfilled with only the orbiter 
engines and aerodynamic surfaces. 

At around 67 seconds the aileron began to deflect. It has not been 
visible up to that time because its cost was high relative to its control 
effectiveness. It was visible from that time on because it was necessary to 
use it to meet the trim constraints, regardless of cost. The limitations on 
the trim power were caused primarily by the rudder hinge moment limits, which 
are shown on the graph. It can be seen that the rudder reached its hinge 
moment bound about three seconds before the control limits were exceeded. At 
that time the other controls began to increase rapidly toward their bounds as 
the rudder was squeezed by its hinge moment limit. The aileron reached its 
hinge moment limit about one second before the infeasible region began. The 
yaw deflection of engine one increased toward its bound and probably reached 
it in the infeasible region. On the other side of the infeasible region (85 
seconds) something similar occurred. The aileron followed its negative hinge 
moment bound for a second and then dipped rapidly to zero. The rudder remained 
exactly on the positive hinge moment limit while engine one yaw deflection 
decreased. 

In the region between 60 and 70 seconds Mach 1 occurs, causing sharp vari- 
ations in the aerodynamic coefficients. In that time span the rudder decreased 
while the yaw displacement of engines two and three increased. This appears 
to be the only consequence of Mach 1 on the control settings. 
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2.3.4 Control Variable Proximity Penalty Functions 

To show the feasibility of using penalty functions to prevent placing 
the control devices on their limits, a penalty function was devised as shown 
below: 




i=l *^6 . 

1 


i=3,5,7,9, . . «n 


Ti[l - 


+ I 


n /26. - 

1 max min 


max min 


(18) 


The last term in the cost function represents a penalty for approaching the 
variable bounds . This was modified in the **COST** and SCOST subroutines 
(Appendix B) and the solution was obtained for the no SRM TVC case at 71 
seconds. The table below compares these results with the pure thrust/drag 
minimum results from Figure 2-7 . 



6a 

6R 

6el 

6e3p 

6e2y 

(Nt) 

Loss 

Thrust /Drag Minimum 

-6.9 

7.5'^ 

3.14 

2,1 

.66 

74000 

Proximity Penalty 

-6,2 

6.7 

4.75 

2,1 

.86 

74500 


6e2p = -6e3p 
6e3y - 6e2y 


^means on boundary (hinge moment limited) 

The penalty function forces the rudder off its limit, but consequently it in- 
creases displacement of the engines (6ely and 6e2y) . It is clear that the 
altered penalty function should cause more thrust /drag loss and the table shows 
that it does. Or course a proximity penalty function does not increase the 
vehicle’s trim power, but it does keep the control devices off their bounds 
until the boundary value(s) is absolutely required for trim. 


2.3.5 Minimum Thrust/Drag Losses With Added Control Power 

Figure 2-8 represents the results of minimizing the thrust and drag losses 
(subject to the mean wind) with SRM TVC in the control loop. As can be seen, 
the aileron did not deflect at all and the rudder only deflected between 60 and 
70 seconds. This seems to contradict the results on the previous graph since 
on that plot the rudder dipped at Mach 1 and the engines increased. However, 
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the costs of the rudder and engines are dissimilar functions. In fact, the 
derivatives of the engine costs are functions of the engine settings so that in 
one region the partial with respect to rudder may be larger than the engine 
partials, and in another region the partials with respect to engines may pre- 
dominate. The solid engines are utilized at around 1-2 degrees deflection, a 
region in which their thrust loss is very small. So simplex selected the rudder 
to use in place of increasing the solid engines and causing a larger thrust loss. 

Figure 2-9 shows the time history of the vehicle with beefed up hydrau- 
lics. (The aileron and rudder hinge moment limits removed) as can be seen, 
the rudder peaked at 12.8 degrees, at 80 seconds. At that time the hinge 
moment limit is actually 5.23 degrees. By allowing the rudder to reach the 
values shown on the graph, the aileron was not required to exceed its nominal 
hinge moment limit. Therefore, the aileron's hydraulic power would not require 
beefing up if the rudder's hydraulic power were unlimited. The character of the 
curve up to 90 seconds is the same as Figure 2—7 , since the hinge moment limits 
do not impact the solution up to that point. 

Figure 2-10 is a plot of the thrust-drag losses versus time for the con- 
figuration with and without TVC and with beefed up aero surface hydraulics to 
give unlimited hinge moments on the aileron and rudder (stop limits remain 
unchanged) . The thrust-drag losses were lowest with TVC and highest with the 
current aero surface hydraulics and no TVC. Of course this plot does not 
reflect the weight increase required to hydraulically control the solid engines. 
Nor does the chart show the weight cost of beefing up the aero surface hydrau- 
lics. As might be expected, the losses were highest when the control power 
approached its limits, on the no TVC case. It is easy to see that when the 
control settings are dictated almost entirely by the trim requirements, little 
can be done to minimize the cost function. Specifically, the high costs 
occurred because the rudder reached its hinge moment limit and the orbiter 
engines (primarily were forced to take up the time slack. As long as 

the orbiter engines were not required to vector more than two or three degrees 
their thrust loss was very small (see Equation (16)). When they were forced 
outside this range their thrust loss became significant. Furthermore the 
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aileron was forced into the trim solution to counteract roll torques and the 
aero surface drag loss was also increased. 

2.3.6 Use of Control Devices to Eliminate Side Force 

It was further desired to assess the ability of the control devices to 
function as force canceling devices. The usual philosophy behind control 
selection and placement is moment cancellation. Almost all of the report up 
to this point has dealt with nulling torque disturbances while minimizing a 
cost function (usually thrust/drag loss). This section imposes the requirement 
of zeroing side force as well as torques. Adding this requirement will theo- 
retically cause the vehicle to fly the exact preprogrammed trajectory regard- 
less of the disturbance. 

It was stated in subsection 2.2 that zeroing side force was probably 
impractical. However, it is impossible to tell the exact control limitations 
using Lagrangian multiplier techniques. Since the Simplex algorithm will give 
a solution if one is feasible, this problem was solved again using the new 
technique. 

Figure 2-11 shows the results of this study, using the mean crosswind 

for the disturbance. As can be seen, the controls met the trim requirements 

until around 65-70 seconds; at that time the control devices exceeded their 

limits. The limiting devices are primarily the rudder and the solid engines. 

When 6 , and 6 approached their 10-degree maximums they caused 6 , and 
24R . e5L 

6 to rapidly increase also. At some time between 65 and 70 seconds the 
e5L 

full control trim capacity was reached. At about that time the rudder also 
reached its hinge moment limit. The limitations on one of these variables 
cause the trim capacity to be exceeded. The aileron and the orblter engines 
do not limit the solution to the extent the rudder and solid engines do. The 
orbiter engines appear to be well within their hard limits up to the time the 
rudder and solid engines reach their boundaries. The aileron is increasing 
quickly but it appears to be well within its hinge moment limit at 66-68 
seconds. 

It is interesting to note that the rudder deflection is negative for this 
case. The rudder deflection was positive when only roll and yaw torques were 
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Figure 2-11. MINIMUM THRUST-DRAG SOLUTION L = N = Y = 0 TVC 
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being cancelled. This happens because the side force requirement dominates 
the solution. The biggest need is for a negative side force, caused by the 
wind, and the rudder furnishes part of that. On the other hand, the aileron is 
positive instead of negative which it was in the previous cases . Apparently 
the aileron, which has no effect on side force, is being used to cancel some 
of the torques created by the deflections of the other devices and by the wind 

torques . 

Judging from the control trim limitations alone, it is evident that using 
the control power to cancel side force is impractical. Other information 
obtained from the study shows that this requirement is impractical from the 
standpoint of thrust and drag losses as well. Figure 2-12 shows the losses 
incurred by trimming the side force are very large compared to those incurred 
from trimming the roll /yaw torques only (Figure 2-10). Note that the curve 
seems to be increasing asymptotically as it nears the infeasible region. This 
is caused by the lack of surplus control authority which limits the flexibility 
of the minimization scheme. As it nears the limits of feasibility this curve 
has the same appearance as the curve on Figure 2-10 which represents no TVC. 

The two curves behave similarly for the same reason. When the trim require- 
ments are very stringent, the minimization program has no room to adjust the 
controls for minimum loss. 

If the aerodynamic hinge moment were increased it would be possible to 
trim the vehicle’s side force. Although only one time point was examined (75 
seconds) it was found that the solution at that point required about 17 degrees 
of negative rudder and five degrees of SRM displacement. The orbiter engines 
and the aileron were negligibly displaced compared to the rudder and solids . 
Experience gained from the study indicates that all other time points will 
trim if 75 seconds has a significant trim capability margin. Nevertheless, the 
thrust-drag loss was around 200,000 newtons, which would seem to be unacceptably 
high. 

2.4 DYNAMIC SURVEY 

In the study to this point the disturbance has been the mean crosswind. 

In order to assess the interaction of the dynamics with the control system 


2-31 



TR'“1246 


NORTNIIOP SERVICES, IRC. 



Figure 2-12 MINIMUM THRUST DRAG LOSS VERSUS TIME ZERO SIDE FORCE MEAN 
CROSSWIND DISTURBANCE 
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command distribution a quick-look dynamic response was generated for the 
NAR configuration. This configuration was chosen because the simulation of 
it was readily available. The data for this configuration are presented in 
Appendix D. Also in this appendix is a block diagram of the dynamics model and 
NAR control scheme. The disturbance was a synthetic crosswind with a gust of 
7,5 m/sec superimposed on a 75 m/sec wind. 

The results of the NAR simulation are shown in Figure 2-13. The commanded 
control settings are the same as the actual control settings if actuator 
dynamics are neglected. As can be seen from Appendix D, the rudder, 6 ^^ and 
5 e 23 D used solely to control roll. Yaw is controlled by 4 ^ 23 l 
solids (d gL 4 *^eL5^ are controlled by feedback from both yaw and roll. The 
torque outputs from these control settings were split into roll and yaw to 
generate a time history of roll/yaw torque commands. 

These torque commands were input to the simplex program. to generate a set 
of minimized trim solutions. The results of the minimum thrust-drag solution 
process are shown in Figure 2-14. The rudder is more active than the NAR 
control scheme allows (no hinge moment limit was known or used) . The top 
orbiter engine ( 6 ^^^ receives a larger displacement and is yawed completely 
to its limit by five seconds. The solid engines and 6 ^ 23 ]) slightly less 
active than the NAR control scheme calls for. There appears to be no signif- 
icant differences in the deflection of 6 . 

The most striking result of this graph is the similarity in results 
between the NAR control system and the Simplex method of torque command dis- 
tribution. It is probable that NAR performed some similar sort of optimization 
to determine their signal-splitting system. 

Figure 2-15 shows the thrust-drag losses for the two control schemes. 

The simplex scheme yields about 10 percent less thrust-drag loss. However, 
this is offset by the fact that the simplex scheme places the top orbiter 
engine on its limit at the time of peak disturbance. 
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The penalty coefficient on orbiter engine one could be adjusted to move 
it off its limit. Alternatively the cost functional could be modified to 
penalize proximity to the control limits. This type of penalty function was 
discussed earlier in the report. 

2.5 FURTHER APPLICATIONS 

There are several possible additional uses for this algorithm. On the 
current NAR configuration there are still payload and load problems during 
ascent. These problems may require changing the control system to resolve 
them. The solution algorithm shbuld give an indication of the way to proceed 
in changing gains, and signal splitting. If soft limits are desired, the 
algorithm could be used to generate position dependent control gains to achieve 
the same purpose. Ultimately it may be possible to implement the system on- 
board. The problem that must be solved here is computer time. It requires 
three to nine matrix row or column operations to get from one flight state 
to the next if the time change is small enough. This system would automatically 
solve an engine malfunction problem, minimizing the deterioration in mission 
capability . 

The algorithm has several uses applicable to new configurations. If a 
configuration does not have the control system defined, a typical torque com- 
mand profile, such as the mean wind, may be used with the algorithm to compute 
the settings of the control devices. This will give a preliminary control 
assessment to determine if controllability problems exist and how to proceed 
with signal splitting. If they are not already set, the algorithm can be used 
to set engine cant angles at their most efficient values. Use of the algorithm 
on problems such as these will add to the background and feel of the designer. 
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Section III 
CONCLUSIONS 


Lagrangian multiplier techniques may be used to solve quadratic cost 
functions only. They also provide a method for blending the commands 
to the control devices. However, the trim problems must be unbounded 
for a solution to always be available. 

The Simplex based algorithm works for general problems .The cost 
function may be of any form as long as it is piecewise differentiable, 
The algorithm considers the variables’ bounds explicitly in the 
solution. 


The algorithm is most effective where there is adequate surplus trim 
power to give flexibility to the solution. 

The algorithm can be used to solve problems such as optimally 
trimming the shuttle ascent vehicle or the orb iter in reentry. 

(’’Optimal" trim means approximately the optimal trim solution wit 
respect to a particular cost function.) It is also useful in obtaining 
a feel for the trim range and an insight into the basis of signal 
splitting schemes on the conventional Shuttle ascent control systems. 

The algorithm may be used in flight if sufficient digital computer 
speed is available. On-board it would solve the engine malfunction 
problem automatically. 

The 049 ascent vehicle has controllability problems if the solid engines 
are not vectorable. If the solid engines are used the controllability 
problem vanishes. Without SRM TVC simplex shows that there is no 
feasible roll, yaw trim solution between 74 and 85 seconds of flight. 
With beefed up rudder hydraulics there is a solution. 


• The dynamic survey performed on the NAR configuration showed the 
simplex control scheme to be close to the results obtained ftom the 
NAR control scheme . The simplex scheme saved about 10 percent of the 
thrust— drag losses compared with the NAR control scheme. 
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ENGINE FORCES AND MOMENTS 


T 

6.68 m 


© 

0 9.34. 


Centerline 
of HO Tank 


eg 


6 Positive 
ep 




^ ey 


6 Positive 


t 

Z 


Engine 1 


F cos 18° 6 


eyn 


force in Y 


(9.34 + Z )F cos 18° 6 roll 

eg ey 


-X F cos 18° 6 
eg ey 


yaw 


where X , Z > lengths 
eg eg 


Engine 2. 


F{cos 12° cos 3.5° 6 + sin 12° sin 3.5° 6 + cos 12° sin 3.5°} 


ey 


ep. 


(6.68 + Z ) (Y force} - 1.345F{sin 12° - cos 12° 6 } 

eg ep^ 


- X F {sin 12° sin 3.5° 6 + cos 12° cos 3.5° 6 }{yaw} 


eg 


ep 
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Engine 3. 

F{cos 12® cos 3.5® 6^^^ - sin 12® sin 3.5® - cos 12® 

(6.68 + Z ) (Y Force) + 1.346F{sin 12® - cos 12° 6 

eg ep3 

- F {-sin 12® sin 3.5® + cos 12® cos 3.5® 


sin 3.5®} (Y force) 
roll 
yaw 


A ventral fin was added to the configuration to improve the aerodynamic prop- 
erties in the region of max q. It significantly reduces the yaw and roll 
moments due to side slip. This effect is accounted for by direct addition 
of AC , AC and AC to their respective terms in the equations. 
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DATA 


FLIGHT 

TIME 

SEC 



NEWT 

THRUST/ 

ENGINE 

ZCG 

m 

■ 

NEWT/m^ 

DYNAMIC 

pressure 

Aref 



bref 

m 

NEWT 

THRUST/ SOL. 

25 

1 .65 X 10® 

-1,58 

23.345 

.482 X 10^ 

3.177 X 

10^ 

2.832 


1 .28 X 10^ 

40 

1.76 X 10® 

-1 .5847 

23,42 

.987 X 10^ 





1.12 X lO' 

50 

1 .825 X 10® 

-1.5914 

23.47 

.134 X 10® 





1 .137 

60 

1.885 X 10® 

-1 .5953 

23.52 

.174 X 10® 





1 .157 

65 

1 .92 X 10® 

-1,5979 

23.545 

.194 X 10® 





1 .166 

70 

1 .94 X 10® 

-1 .60 

23.57 

5 

.212 X lO'’ 





1.175 

75 

1 .97 X Id® 

-1 .4626 

24.13 

.226 X 10® 





1.183 

80 

1.98 X 10® 

-1 .455 

24.18 

.233 X 10® 





1 .192 

90 

2.025 X 10® 

-1 .440 

24.33 

.217 X 10® 





1 .205 

100 

2.04 X 10® 

-1 .4327, 

24.535 

.165 X 10® 





1 .213 

no 

2.06 X 10® 

-1.4255 

24.74 

.117 X 10® 




f . 

1 .216 

145 

2.07 X 10® 

-1 .400 

25.62 

.231 X 10^ 


f 



1 .0 


1-3 









00 FT^ - FIN 



DATA IS PER/DEG 


FLIGHT 

TIME 

AFT 

ACmB 

■ 

^1. 

FORWARD 

ACmB 

ACyB 

25 

,0064 

.0031 

^.004 

-.011 

40 

1 

,0067 

.0032 

-.0044 

-.012 

50 

,0074 

.0033 

-.0048 

-.013 

60 

.0085 

.0036 

-.0056 

-.015 

65 

.0094 

.0038 

-.006 

-.016 

70 

.014 

.0042 

-.006 

-.017 

75 

.01 

.0 

CM 

o 

-.0058 

-.0165 

80 

.0088 

.0035 

-.005 

-.014 

90 

.0075 

.0027 

-.0044 

-.0105 

100 

.005 

.0017 

-.0028 

-.008 

no 

.004 

.0014 

-.0022 

-.006 

145 

.0028 

.001 

-.0015 

-.004 
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DATA 



^9 

TOTAL 

VELOCITY 

m/sec 

VR 

Uo 

m/sec 

liJO 

m/sec 

Cos 

S1n 

GRAVITY 

g 2 

m/sec 

MACH § 

WIND 

VELOCITY 

m/sec 

■ 

1 .2*= 

25 

.954 X 10^ 

.954 X 10^ 

.279 

-.012 

1.0 

9.8 

.28 

2.0 

.0007 

3.4® 

40 

.150 X 10^ 

.149 X 10^ 

.165 X 10^ 

.050 

.999 

9.8 

.44 

9.0 


4.5 (4.6) 

50 

.190 X 10^ 

.186 X 10^ 

.340 X 10^ 

.125 

.992 

9.79 

.56 

15.0 


5.7 (5.9) 

60 

.241 X 10^ 

.232 X 10^ 

.654 X 10^ 

.222 

.975 

9.78 

.707 

24.0 


6,1 (6.5) 

65 

.272 X 10^ 

.257 X 10^ 

.875 X 10^ 

.274 

.962 


.8 

29,0 


6.4 (6.9) 

70 

.305 X 10^ 

.283 X 10^ 

.115 X 10^ 

.329 

.944 


.9 

34.0 


6.7 

75 

.343 X 10^ 

.310 X 10^ 

.148 X 10^ 

.384 

.923 

9.77 

1 .01 

40.0 


.6.5 

80 

.385 X 10^ 

.337 X 10^ 

.187 X 10^ 

.449 

.893 


1.13 

44.0 


3.53 

90 

.486 X 10^ 

.392 X 10^ 

.288 X 10^ 

.566 

.824 

9.76 

1.43 

30.0 

.007 

0 

TOO 

.612 X 10^ 

.445 X 10^ 

.420 X 10^ 

.664 

.748 

9.74 

1.8 

0 



110 

.768 X 10^ 

.498 X 10^ 

.584 X 10^ 

.681 

.732 

9.73 





140 

.152 X lO'^ 

.673 X 10^ 

.126 X 10^ 

.874 

.486 

9.68 





I 


N> 

O' 
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DATA 



“■Sr 

■''Sr 

C„ 

"sr 




q 0^ 
^0 0 

25 

.273 

.504 

-.510 

-.283 

.302 

-1.66 ^ 

-.203 X lo’ 

40 

.265 

.408 

-.489 

-.285 

.315 

-1 .68 

-.512 X 10"^ 

50 

.259 

.462 

-.473 

-.286 

.325 

-1 .70 

.450 X 10'^ 

60 

.215 

.394 

-.388 

-.291 

,404 

-1 .83 

.112 X 10"^ 

65 

.181 

.319 

-.310 

-.298 

.468 

-1 .99 

-.158 X 10*’ 

70 

.173 

.300 

-.345 

- .326 

.460 

-2.05 

.687 X 10“^ 

75 

.206 

.292 

-.340 

-.384 

.344 

-1.97 

-.103 

80 

.186 

.217 

-.254 

-.356 

.266 

-1.92 

.412 X 10’^ 

90 

.105 

.132 

-.137 

-.299 

.238 

-1.93 

-.277 X 10"^ 

100 

.055 

.961 X 10'^ 

-.105 

-.246 

.269 

-2.03 

.301 X 10'^ 

no 

.0406 

.749 X 10'^ 

-.077 

-.196 

,207 

-1.98 

-.851 X 10'^ 

145 

.0286 

.573 X 10'^ 

-.061 

-.122 

-.0284 

-1 .60 

-.100 X 10'^ 


-6 
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DATA 


lilB™ 

c 

■"sa 

'Sa 


25 

.0458 

-.0430 

.218 X 10 ^ 

40 

,0444 

-.0458 

.201 X 10 ^ 

50 

.043 

-.0487 

! 

.190 X 10 ^ 

, 60 

.0358 

-.0544 

.179 X 10 ^ 

65 ^ 

.0344 

-.0630 

.174 X 10 ^ 

70 

.0301 

-.0630 

.169 X 10 ^ 

75 

.0258 

-.0544 

.160 X 10 ^ 

80 

.0244 

-.0458 

.154 X 10 ^ 

90 

.0172 

-.0286 

.144 X 10 ^ 

100 

.00286 

-.0215 

.133 X 10 ^ 

no 

-.00286 

-.0158 

.122 X 10 ^ 

145 

-.0114 

-.00859 

.914 X 10 ® 


TR-1246 


Ix „ 
Kg - 

Iz 2 
Kg - 

.953 X 10 ® 

.591 X 10 ® 

.856 X 10 ® 

.547 X 10 ® 

.794 X 10 ® 

.519 X 10 ® 

.733 X 10 ® 

.491 X 10 ® 

.702 X 10 ® 

.478 X 10 ® 

.671 X 10 ® 

.464 X 10 ® 

.629 X 10 ® 

.420 X 10 ® 

.606 X 10 ® 

.405 X 10 ® 

.559 X 10 ® 

.375 X 10 ® 

.512 X 10 ® 

.346 X 10 ® 

.466 X 10 ® 

.317 X 10 ® 

.329 X 10 ® 

.228 X 10 ® 









NORTHROP SERVICES, INC. 


TR-1246 


DATA 


FLIGHT 

TIME 

RUDDER HINGE 
MOMENT LIMIT 
DEGREES 

<^Sa (INCLUDED ANGLE) 1 

^ AILERON HINGE 
MOMENT LIMIT 
DEGREES 

25 

No Hinge Moment Limit (30) 

No Aileron Limit 

40 

42 (30) 

71 .8 

50 

30.8 (30) 

52.6 

60 

23.5 

40.0 

65 

14.7 

25.1 

70 

8.19 

14.1 

75 

5.54 

9.47 

80 

5.23 

8.91 

90 

6.27 

10.69 

100 

10.23 

17.5 

no 

19.67 

33.64 

140 

No Hinge Limit 

No Hinge Limit 


Hard Limits 

Hard Limits 


(±30") 

(40° up - 15° down) 


The hinge moment limits are a function of dynamic pressure and as such 
would change with trajectory variation. However, this data assumes a nominal 
ascent. In addition, I have included the hard deflection limits which govern 
when the hinge moments do not. I do not have any data for it yet, but there 
should be a yaw introduced by unequal aileron deflection. 
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The simplex algorithm is a well known algorithm for solving linear, 
bounded problems with linear cost functionals. It is uniquely suited to the 
trim optimization problem since simplex solves the problem of more variables 
than constraints and also handles explicit inequality constraints. The orig- 
inal simplex algorithm is as follows: 


1. Examine the a., j = 1, n-m. 

J 

a. If all a° > 0, the solution is optimal. 

b. If some a? .< 0, choose j = s by a° = min a° < 0. 

J s J 

1 

2. Examine the a,l=l,...,m. 

s 

a. If all ^ 0, the value of Z is unbounded in the class of 
solutions given by (7). 

b. If some > 0, choose i = r by 


i i 

min a >0 — 

s a 

s 


^ >0 
r 
a 

s 


(b^ > 0) 


3. Using a^ for a pivot, perform one elimination to obtain the new 

canonicll form. The elimination is applied to the right hand side 
to give the new values of the new basic variables and an extra E.R.T, 
applied to row zero to give the new relative cost factors. The new 
value of Z (with sign changed) appears in row zero of the right side, 


4. Return to 1, 


B-1 


( 
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It can be shown (ref. 2) that if a feasible solution exists, we can obtain 
a basic feasible solution. Applying the simplex algorithm to this solution, 
assuming Z is bounded below, we either prove the solution optimal or obtain an 
improved basic feasible solution, provided all b^ > 0 for > 0. Thus we can 
never return to a previous solution since each one has a lower value of Z than 
the preceding one. But there are only a finite number of basic solutions (not 
exceeding the number of combinations on n things taken m at a time) so that we 
must eventually arrive at some best basic feasible solution. But this must be 
an optimal solution or we could repeat the algorithm. Hence there exists a 
basic feasible solution which is optimal. This proves that if any feasible 
solution exists and Z has a finite minimum then there exists a basic optimal 
solution. Furthermore it shows that the algorithm is a finite iterative pro- 
cedure, i.e., we arrive at a solution in a finite number of steps. This is 
also clearly true if Z is unbounded below. 

Finding A Basic Feasible Solution - The number of linearly independent 
rows of a restraint matrix is called the rank of the system. If the rank p is 
less than m, then either there are m-p redundant equations or else there is no 
solution at all, feasible or otherwise. If the rank is m, then there is some 
basic solution but there may be no feasible solution. Hence, given a problem 
there is, in general, no way of knowing whether any solution exists. 

To overcome these difficulties, a very simple device exists which either 
finds a feasible solution or determines that none exists. It also induces 
rank m in case p < m. The idea is to embed the given problem in a larger 
system and solve a preliminary LP problem before attempting to minimize Z, 

This process is called Phase I of the simplex method, the optimization of Z 
being referred to as Phase II. 

We suppose that all b^ ^ 0. (If not, the entire equation for b*" < 0 can 
be multiplied by -1.) 

We now define m artificial variables 


n+1 

x 


> 0 , 


n+2 

X 


> 0 , 


n-fm 

X > 


0 


B-2 
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and which will eventually be restricted to be non-negative. The co- 

efficients of these variables will be new vectors arbitrarily adjoined to the 
restraint matrix. The coefficient of will be unity in row i and in an 

added row w, zero elsewhere. Letting m + 1 = w the variable x*^ will have a 
coefficient only in row w. The system of restraints is now 

I x"-"^ = 0 


i=l 

y a^ x^ + x'^^^ = b^ (i = 1, . . . , m) 


(B-2) 


Since the x^^^ are artificial, they must not appear in the final basic solu- 
tion, except perhaps at zero level. If > 0, then (B-1) guarantees that 

they are all zero since a sum of non-negative numbers equal to zero implies 
each is zero. 

If the = 0, then clearly (B-2) is the same as the original restraints. 

Note that the rank of the system (B-1), (B-2) is m+1. This is clearly seen by 
displaying the entire expanded tableau: 


1 

X 

2 

X . . 

n 

. X 

n+1 

X 

n+2 

X 

n4m 

... X 

n+w 

X 


_ 

o 

^1 

o 

^2 * ’ 

o 

. a 

n 

0 

0 

. . . 0 



Z 

1 

1 

^2 • * 

1 

. a 

n 

1 


... 0 




• 

• 

* 

• 

• 




• 

* 

• 


• 

• 




• 

• 

• 

♦ 

• 


• 


— 

* 

* 

* 

• 

• 


• 



• 

* 

• 

• 

• 


• 



• 

m 

^1 

m 

^2 * * 

m 

. a 

n 

0 

0 

... 1 



b“ 

0 

0 . . 

. 0 

1 

1 

... 1 

1 


0 
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Although (B-3) is not quite in canonical form, there is an obvious starting 
solution if we allow to be negative, namely: 

x^=0, j=l, — ,n 


n+w 


m . 

= -I 


i=l 


n+i , i . T 

X =b,i=l, 


(B-4) 


The solution (albeit infeasible) will remain valid of we apply the following 
E,R.T: 

a. = a. - a. (j=l n+m) , b = b - b 

J J J 

i.e., if we subtract Equation 1 from Equation w. We can then subtract equa- 

w w 

tions 2, 3, ...» m in turn from Equations w. If we redefine a^ and b as 


w 


m 


i = I (j = 1, . ... n+m) 


j iSl J 


m 


b” = - I 

i=l 


the result is 

1 

X 


the following tableau: 
2 n n+1 

X ... X X 


n+2 

X 


n4w 

X 



(B-5) 
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w H"tw < 

which is in canonical form except that b <0, and hence x <0. The first 
objective will therefore be to maximize x”^. Zero is clearly an upper bound 
for x*^^ by (13). On the other hand, any feasible solution to the original 
problem would leave all x =0. Therefore if max x <0, there is no 
solution to the given problem. 


From Equation w in (B-5) , it is clear that maximizing x^~^ is the same as 

n . 

r w i 

minimizing F = I a. x since their sum is a constant, 

j=i ^ 


n-fw 

X 


+ F 


= b 


w 


w 

Hence we simply use the a . for 

o ^ 

to be ready for Phase II if 


relative cost factors in Phase I, carrying along 
and when x^^ = 0. 


During Phase I, maximizing x , its column of coefficients must never be 
a candidate for elimination. 


It is not possible to eliminate all the artificial variables from any 
basic solution even in Phase II. For even if the m original restraint rows 
a^ (i = 1, ...» m; j = 1, . , . , m) were linearly independent, row w is a linear 
combination of them and hence the m + 1 rows for i = 1, ..., w and j = 1, ..., n 
are linearly dependent. Hence at least one artificial column must always remain 
and we have insisted that column n + w stay in during Phase I. If, however, 
the original m rows had rank p < m, then m - p + 1 artificial columns will 
remain at the end of Phase I. If x^"^ =0, the solution is a valid one; if 
^n+w < there is no solution. 

niw 

Degeneracy may also cause x to reach zero before the other artificial 

columns are removed from the basis. In this case, some of the aY may still be 

n4w 

negative. Thus Phase I should be terminated when either x = 0 or when all 
w 

a. 0, whichever happens first. 
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It is possible thet x will be eliininated in Phase II if some otheir 
artificial column is still in the basis. Whenever any artificial column is 
eliminated from the basis, it is simply dropped from further consideration. 

If the original restraints contained some legitimate unit vectors, then 
the corresponding artificial columns need not be added and the corresponding 
rows are not subtracted from row w. 


B-6 
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This appendix Is a listing and description of the computer program 
developed to solve the problems in this study. Figure C-1 is a flowchart of 
the calling sequence for program SMPLX. 

DATRD is the program that reads the data and calculates the control 
devices* stability derivatives. This program also reads in such things as 
penalty coefficients and initial solution and starting range. The program 
can be started with a zero vector for the starting solution and 100 percent 
range to insure a feasible solution if one is available. 

Then the program enters SIMTM which calls all the subroutines that 

calculate the solution. For the first pass the partial costs are set 

9 Z 

(subroutine cost calculates at the initial solution, the temporary 
bounds are set up about that solution and the simplex tableau is set up and 
solved. This gives the minimum cost sum (subroutine SCOST) on the temporary 
range if the costs are linear enough on that range. 

On subsequent passes CVARB calculates the actual variables from the 
simplex transformed variables, TORQ calculates the sum of the variables 
times their coefficients and the disturbance, and SCOST calculates the 
actual total cost so that actual Improvement can be discerned. Then the new 
cost is compared to the old cost and if there is no improvement the bounds 
are squeezed by half (SQZ) . If there is Improvement the new partial costs 
(SCOST) are calculated at the new solution. Then the new temporary bounds 
are set. Subroutine SETUP puts the problem in simplex form and SOLTN solves 
it. The improvement criteria is then checked and if it is fulfilled, the 
program loops again. 

The only things that require change to solve a different problem are 
COST and SCOST, and then only if the cost function changes form from Equation 
(16). SCOST represents the total cost and if the form is changed, then the 
cards from "TCOST = 0.0" to "RETURN" must be altered according to the actual 
cost equations. Then in subroutine COST the partial derivatives of the total 
cost equation, evaluated at the current solution, are calculated. To change 
COST, simply remove the cards from "18 CONTINUE" to "DO 22 I = 1, NVAR". 


C-1 
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Fiqure C-1 . CALLING SEQUENCE FLOWCHART OF PROGRAM SIMPLEX 
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The number of variables in the problem is represented by Fortran mnemonic 
"NVAR". The number of equality constraints is called "NEQ". These must be 
changed if the problem size changes. The program is dimensioned for a maximum 
of three equations and 12 variables. 

An input description is listed below: 


CARD 

FORMAT (5E13. 4) 


1 

YEl 

ZEl 

PAEl 

YAEl 

2 

YE2 

ZE2 

PAE2 

YAE2 

3 

YE 3 

ZE3 

PAE3 

YAE3 

4 

YE4 

ZE4 

PAE4 

YAE4 

5 

YE5 

ZE5 

PAE5 

YAE5 


(YEl is y position relative to eg of engine 1, ZE2 is z position engine 
2, PAE3 is pitch angle engine 3, YAE4 is yaw angle engine 4) 

FORMAT 13 

6 TIME (time of flight) 

FORMAT (5E13.0) 

7 THR, ZCG, XCG, AREF, BREF 

2 

(orbiter thrust per engine (nt) , ZCG 9m), XCG (m) , reference area (m ), 
reference length (m)) 

8 QBAR, 6, MASS, IX, IZ 

2 2 

(dynamic pressure (nt/m ), gravity (m/sec ) mass (kg), x axis inertia 
(kg m^) , z axis inertia (kg m2)). 

9 THRS , BELT 

(solid engine thrust per engine (nt) , initial solution range, rates 

of full range 0-1.0) 

10 (DIMAX(I), DIMIN(I), I=1,NVAR) 

(maxima and minima for each variable one at a time) 

11 DCCB, DCNB, DCVB 
(A AC , AC ) 

*■6 ”b ^6 

12 CLB, CLDA, CLDR 
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CARD 

13 


14 


15 

16 

17 

18 


CNB, CNDA, CNDR 

(C , C . C ) 

n ’ n n. 

CYB, CYDA, CYDR 

(C , , C ) 

^ 'a 'r 
(TC(I), 1=1, NEQ) 

(disturbance torques and forces aside from wind) 
BETA 


(wind velocity divided by vehicle total velocity) 
(PC(D), I=1,NVAR) 

(penalty coefficients) 

(A(I,KKK), I=1,NVAR) 

(initial solution vector (temporary location)) 


C-4 



// «SHARD X3AR 0250 


// FOR 

■jjiONF WORD IMTFG 


*mCS U 13?DRINTER#CAR0 ) 

*OEl t*\PASS) 


INTEGFP TI^’T 
REAL MASS i-IX.IZ 


COMMON A(17*?6 ) »W{ 17 ) tLClT) *I 1 T*JJ*I I *IW . 

MMOM DIMAX( 12 ) *DIMIiM ( 12 } *PC(1?) .00(12) 


COMMOf^ FIMAX(12) .FIMIN(12) *C(3*13) .TCO) . DFLT . TCOST ♦ COST 1 . COS r 2 

7,2957795/ .INVAR 


COSTl=10O0Q00, 
0 



READ (NR.200) YE 1 , ZE 1 . PAE 1 » YAE 1 
READ (MR. 20 


READ (MR. 200) YE 3 . ZE3 . PAE 3 ♦ YAE3 
(MR. 2 


READ (MR. 200) YE 5 . ZE5 . PAE 5 ♦ YAE5 
ni TIM 


READ (MR. 201 ) ThR .ZCG . XCG . AREF .BREF 
. 201) QBAR.G.M 


READING, 201) THRS.DELT 

( T ) . 


READ (MR. 201) DCLB.DCNB.DCYB 
( MD » 2m 1 ri R .ri I 


READ(MR.20l) CNR . CNOA . CNDR 


rasK 




I ri 9i.i M 


lT4 Sn«l> 4 ft 7 iTMriH: 


REA0(MR.201 ) ( TC ( I) ♦ I = 1 . M EQ ) 
) ! 


FZCG s ( ZCG+1.47 ) /BREF 




= QRAP*AREF 

* 



ENG I ME 1 PITCH DOES MOT CONTRIBUTE TO LATERAL FORCES AMD MOMENTS 


TEMP s PAFl/CMV 

= COS (TEMP) 


= PAF2/CMV 




SPE2 

TEMP 

SIM (TEMP) 
= PAE3/CMV 

CPE3 

SPE3 

= COS (TEMP) 
= 5 IN (TEMP) 

TEMP 

CYE2 

' = YAF2/CNV , 

= ros(TEMP) 

SYE2 

TEMP 

= SIN (TEMP) 
= YA=^3/CNV 

CYF3 = COS{TEMP) 
5YE3 • = SIN(TEMP) 


YAEA/CMV 
s ( TFMP ) 



TEMP = YA^5/CMV 

CY 


CLEIY = (ZEI+ZCCO^^THR-WCPEI/OSB 

1 .n*xCA*THR*CP.El/Q5B 


CLE2° = -(YE?*CPE2 + (2CG + ZE2)^^SPE?*SYE2)*THR/QSB 

rNF?P ° 

CME2P = -1 .0<^XCG^^THR*CPE2^^CYE2/QSB 

) »THP*rPF2»CYE2/QSB 


CNE2Y = -1.0*XCG<^THR*CPE2^^CYE2/0SB 

»TMR»rPE3*CYE3/QSB 





CLE3Y « 
rWF3Y = 


CLEAR s 
CLFAL « 


CMEAR = 


CLE5L * 
MEAL = 


CNEAR = 
CNEA 
CME5P = 
CN 

CNE5L =' 


(ZE3+ZCG )*THR-»^CPE3*CYE3/QSB 
-N- 


0«707*THRS*(-YEA- { ZEA +ZCG)*CYEA) /OSH 
0>707»THRS» (-YEA-MZEA 4-Z CG ) »CYE A ) /QSB 
-0«707*XCG*THRS*CYEA/OSB 
r),7r>7»THRS»(^YE5-UE5 ->-ZCO ) 1 !^Y£ 5 ) /QSB 


0.707^^THRS*<-YE5+(ZE5 +ZCG ) *CYE5 ) /OS9 


0. 707*XCG^^CYEA*THRS/OSB 
'0^707*XCG->fCYEA*THRS/QSR 
-G.707*XCG*THRS*CYE5/QSB 
0>707»XCG»CYE5»THRS/QSB 
•0.707»XC6»CYE5-“»THRS/Q5B 
MF A 



COMMUTES MOMENT COEFFICIENT RELATIVE TO FOMENT REFERENCE POINT 



C( 1*6) = CLE2P 
M *7 i = 



CU*10)* CLEAR 
1 2 Ci^ 


C(l*12)= CLFSR 
(1*13)= CLE^ 



C(2*2) « CNDA 

5 « r\tnp >■ CYDR»FXCG 



C ( 2 * 6 ) = 


C(2*e) = CNE3P 










C(?»1()) = CNF4P 
C(?fll )=CMF^L 
^(2tl?)=CNE5R 
r (2 ♦! 3 ) =C''iF/-'L 


C(3tl) s n.o 

£J3f?) = 0.0 

C ( 3 * 3 ) = 0,0 
C (3WQ ^ C’-^ EIP 
C{3*6) = 

^L3.±6 ) ^ CVF?P 
C { ^ ♦ 7 ) ^ 0.0 
C(3»P ) = r‘^F3P 
C ( 3 *9 ) = 0.0 
r(3»10)=C*-^FAP 
C<3*11)= ^^^E■4L 
C(3*12)=CV'F3R 
C(^*i3)= r-vr'=>L 


00 A < = 1 ,VFO : 

TC(K) = TCC^) - C 1 ) 5^-PFTA 

^^■CO^UI^JUE ^ 

REA0(.NP,20 1 ) ( PC ( I ) * 1 = 1 *>:VAR ) 

read (Np» 20 1) (A( I »25) «I=1» WAR) ^ 

wpiTFC "’Vy*P0O) TIME 

WPITE(->W»203) ^1^•AX 

WR I TE ( 'IW ♦?03 ) DIM IN 

WP ! TF ( M'<» .99 ) ^ 

'*'RITF(NW*501) 

WRITE t >^lW . 20 3 ) { C ( 1 ) < I = 1 ♦ a ) ^ ^ ^ 

WRI^F (^'W*'?22) 

WR.I T£ ( A|'>Li203_j ( C ( 1 » n > I =9 ^ NVARl ) 

WR I TF { ^!W ♦ S02 ) 

WRlTF{vW,703)(C(2>I)rI_=l.B) .. 

'/."^ITF f MW* 523) 

TMRITE (MVS?03HC(2*I)*I=9*:'.VAR1) 

WP. I TF ( '■”■'*5 03 ) 

vRITFC-”-'*203)(C(3*I) *1 = 1*8) ^ 

WPITF (MW, 3 24) 

■MRI'f'E (NW*?03) (C(3_*Ti_:._I=9»NVARl) 

WPITF ( M'---'*2 02) (PC( I) »I=1*RVAR ) 

‘-‘'R I T^E ( -w .300) ( TC ( n * I = 1 * NEO ) ^ 

CALL L I ( S IMTV; ) 

99 »^QRA:AIJ 2X.//»54H LATERAL NQN-D IP £NS I ONA L ST ARU. I T V DLi-H V AT I V E.S ^1/ 
rRA0^5 ) 

lO0_'^a_RVAT( 13) 

?on ‘■'ORVAT ( FE13.4) 

201- FORVAT(5F13.0) ^ 

?02 format (BF l 3.3 ) 

inl F-QRM^T(BE13^3) L 

5-OQ FORMAT ( IHl */*10X *6H TIME = *I3*5H SEC $///,) 

501 FORMATt// *5X»5H CLB *8X*5H CLOA*3X*5H ClDR*3X* 

16H CLEXP*7X*6M CLE1Y*7X*6H CLE2P»7X*6H CLE2Y*7X*6H CLE3P) 

502 FCRMAT(// *5X»5H CNB *8X*5H CNDA»8X*5H CMDR*3X* 

16H CN’F1P*TX*6H CNE1Y*7X*6H C-ME2P*7X*6H CNF 2Y * 7X . 6H GME 3P ) 

503 FORMATi// * 5X ♦ 5HCMALP * 8X » 5H CMDA*8X*5H CMDR*3X* 










// J0«=». T 


Lor; drive CART 5^EC CART AVAIL PHY DRIVE 


OOOC 


DOC? 


0002 


0000 


••'? v-10 /^CTUAL ‘"K CONE 10 B 


ff *5(_jAPP X36B ^?B0 


// FOP 

* I OCS ( 1 1B?PRINTFR«CAR.DJ 


WORD INTEGERS 

* A P I T H V E T I r TRACE 



data NW /3 / *NV/AR/ 12/ 

DATA r\t\//B7*2957B/ 


■NVARI =?*NVAP +1 



CALL COST! DO tDIMAX <D I MI N *A*PC) 


D( I ) -A( I *^■VAR1 ) 

• ^ / T ^ = n M 1 


B2 CONTINUE 


, n Pi p r^ p 


COSTl =10000000000000* 

CO TO .12 

CONTINU‘D 

c ^ \7 CVARB ( L *NVAR >A * Df E IMIN ♦ E I VAX *NVAR 1 ) 

TtYAWT*PlTTJ_ 


SCOST ( COST 2 » COSTl •TCOST *D *PC # CN V *D I MAX * 0 1 M I * 


F(COST1-TCOS 


CONT I NUF 
WRITE(NW»2C)?)_ 

MN=1 

TCOST s COSTl 


CALL SOZ . 
CO TO 


CONTINUE 
DO 87 1=1 »NVAR 


DD( I ) =DI I) 




TCOST 


tDELT*TC05T*C0STl ) 


■9 







PA*; 




IOCS 


R^OUIR^VENTS FOR- 

COiM'-'nN' 117 f, VARIABLES - 40 PRvOGRAM 3 22 

1 

f 


E\o 0=^ co ‘- dilation 



C-10 









// 'MJD 

*ST09FCI WS UA SIA'TM 

^^DUf^o IIA rn SliXTM 

// Pnp 

*0~-tF WORD INTPriPRP 

S'JBPOUT SOZ ( DfiLT *TCOST ) 

A = ARS( ^ 

K=AP5( TCnsT ) 

A= ( A~n ) /A 

A = A ^ .S ( A ) 

IF(A-0>01) 10»lv^*15 

15 IF(A-DFLT) 10*10*20 

10 CONTINUE 

DELT=0.5*DFLT 

_.20 CONTINUE 

RETURM 


cr 



A ~ n r 


/ ^ 

.^)Tr,oF 'a'<; ua sqz 

/ ^'09 

nvr WORF) INTFGFRS 

SUF^ROUT I \>F SETUP ( E I MAX « E I r/ I i\ ♦ NVAR t C » TC » A » W t L f RVAR 1 » R'V » I I « J J * H Q 

FUSION ^I^^AX(12) *EIMIN(.12) *C(3*13) »TC(3) *A( 17*?.6) * 'V ( 1 7 ) • L ( 1 7 ) 

pr^EMSIO^t ^HS{3) ^ 

RH5<1) =0.0 

R!-i$(2) =0.0 ■ ^ — 

RH5(3)=0.0 


11 = 18 — — ^ — 

JJ=?3 

nrj A? T = 1 * ^v/AR 


J=I+1 

PHS (11 = RHS ( 1 ) 

-C( 1 .J) *EIMIN ( I) 


A P 

RhiS ( 3 ) =RHS ( 3 ) -C ( ^ ♦ J ) -»r I VI N ( I ) 
oi IS f 3 ) = PHS ( ) -C ( ? t J ) I M 1 N ( I) 


RMS ( 1 ) = Tr ( 1 ) +RHS ( 1 ) 

P H A { ? ) = T r { 7 ) + R H S ( ? ) - — 

RMS { 3 ) =TC ( 3 ) +PHS ( 3 ) 

= RHS ( 1 ) / APS ( RHS ( 1 ) ) 

SI ON 3 = RHS (?) /APS ( RF4S( 2 ) ) 

c.iCMa^=pH,S{ 3 ) /ABS ( RHS( 3 ) ) _ ^ — 


A ( 1 *.NVAR1) = n. 
I I I = I I+l 

0 

13LPROG 


DO 1010 1 = 1 ♦ I I I 
*'.( ( T ) = 0 - 0 


14LPROG 

ISLPROG 

imp 

in )=0 

rOMT T MllF 


16LPR0G 

DO lOOR K=2fIII 

nn inno ^ ^ ^ ^ 

1 noQ 

A(l^ ♦J) =0.0 
rriMTTMitr 



DO 47 1=1 ♦NVAR 

esT-*-! ^ 



C(l#!^ )*SIGN1 

C ( 2 • K ] »S I GM2 ^ ^ 


A 7 

A (4* I ) =C (3 ,<)*SIGN3 

rnMTTA'tiF 


A(2*NVAR1)= RHS(l)i^SlGNl 

A ( 3 . M\/ AP n = RHS ( ? ) S I 2 

A (4 tMVARl ) =RHS ( 3 ) *SIGM3 


RELATION OF 
ANin rOMIROL 

TRANSFORMED VARIABLES ♦SLACKS 
L 1 MIXS ^ 


no 4R 'Csl.'^'fVAR .... 

v'IsK+A 
^J= i^ + NVAP 

UKl) <J 

A { K I . < } =3.0 

4R 

A {< I =1 .0 
DO 46 1=1 .NVAR 



46 

\j=I+4 

j.NVARl) = FIMAXf r )-EI>OlNlI) 


l'P04 

WR I TE ( NW ♦ 1 004 ) 
FORMAT (12H DATA 

. LOADED) . - 

IILPROG 


C-12 







C-13 



C-14 



5 


UA *' CV.^RR 


r 1,-AprN plTEOrr?s 

tTM^FTTC TRACE 


U^POUT I COST ( D * D AX *0 K‘ I M * A * PC) , 

^T^-'TPSTQM nn ? ) ,OTyAX{ 12) 12 ) *A.( W >g6) ♦PC ( 

OAT A NVAR/12/ 

OAT A C^!V/S7« 2937PV — — 


I a=7*MVAR 
00 IP I = 1PIVAR 
J=I+MVAR 
A(l»J)=0>n 


CO NT INUE 
OO IT=1>2 


A(1,IT) = pcnTP^O(IT)/ABS(D(lT) ) 

-1 


NVArv^NVAR 
DO 20 IT = 3P’VAV*2 
TT1rIT+l 


5 _ 50QT (O { IT ) ^<-D { IT ) +0 n T 1 ) ( I T1 ) ) /CNV 

A ^ AIN(S )/(S^CN^Al 


A ( 1 » I T 1 ) - { I T 1 ) -^ P { I T 1 ) * 5 

AM 4 ! T ) =oc ( 1 T ) ( P ) 


DC ?2 T=X*NVAR 
I+NVAR 


A ( 1 ) ss-A{ 1 » I) 






PAr.F on 


// DU^ 

♦ STOR‘D UA *' SCO ST 

// FOP ' 

0 N E W n P I 'V T E f"' E S 

SU‘^~^POUTI^'P TOPO(i^VAR<C»0*TC»ROLLT » 

D I E N S 1 0 M C ( S f 1 3 ) ♦^(12) * T C ( 3 > 

POLI-TgQ. 0 

YAWTsO* 0 

i:--'ITT = 0.0 - 

DO 73 I=1*NVAR 

J=H-1 — 

P0LLT = P0ULT+C{ 1 * J )*E>( I ) 

PITT = PITT + r(3»J I ) 

YAWT = yaV'T + C ( ^ ♦ J ) ( I ) 

9 01 1. T = POl-t T-TC ( 1 ) 

YA'vT = yawT-TC( 2 ) 

D I TT=p \ TT-TC( 3 ) 


73 


> 



•s'ORO I’tTEGfRS 

SUB ! ^ niJTTMF BT[^M(EIMAX^EirMiU NVAR>Dn*nf£i T ♦ D I U-' N ) 

D I MF ^! S ! ON E I MAX ( 1 2 ) • F I iM I N { 1 2 ) » •"»D ( 1 2 ) » 0 I !^'AX ( 1 2 ) * ; I I !'• I 1 2 

00 5 7 ! = lfN!VAR 


EIMAX { I ) 
g- 1 V I N ( I) 


57 


CONTINUE 


I 0^(1 )+DELT*( DINAX( I )-DIMIN ( I )T 
= r)!)( I ) ^DFLT» ( D I MAX ( n-DIMI M ( ! 11 


c 

DO 51 I = 
TFMP = 

ItMVAR 

FT MAX { T ) ->D IMAX (I ) 




IF (TEMt5) 5A »5° >5A 



A 

C T A y ( t ) 

s DTMAXTI) 




) H 

5P 


E I M I N ( 1 ) ~D I M I A' ( I) 




T F f TFVT^) 50 ,51 , 5 1 



S P 

F I V T M ( I ) 
r n T T N 1 1 F 

- n I N/ 1 N ( I ) 




DO 52 1= 

1= T -1 

4 ♦ ♦ 2 




A s S 0 P T { ^ I A X ( I ) * F I ‘■■-'A X ( I ) + E I A X { J ) F I ’ A > ( J ) ) 




T F { A-DT V 

'AyM^)5fi»5(S*53 — — 

— 


5? 

4 p 1 : i ft 1 ~ f 

CONTINUE 




F I .'-'.A X { I) 

FTVAX f T 1*DTMAX ( I) /A 

— 



PI^AXt j) = FIVAX ( J) *0K4AX ( J) /A 

A =: ^q-?T(FIvIN( n^FrUNU )+EI^MN(J)^Er MM(.^ 

s = AB.<( n I M IN ( I) ) 


55 

I t [ A-^ — 

F I T N ( I ) =E I "IN ( I ) *B/A 

CT^/TM< n =F J)*B/A 


52 

= SORT ( F I '^'I -M ( J ) *F N ( J ) -fE I MAX ( I ) *F ( I i > 

T F f A— R i 5 0 • 60 • 6 1 


5 i 

F I M I M { J ) = r IMI N ( J ) / A 

ftmAXCT) = FT''^AX(I)^<'B/A — 


5 0 

A = .SQPT('^I‘’'^i:’(I )^^FI^'lNn )+EIVAX{J)^r,I^^AX(J} ) 

TFfA-Rl 62»62«£^ ^ 

— 


c I 1 1\[ ( I ) = I I N ( I ) ■“ R / A 

cTvAy(n = ri MA X I J ) »B/ A. 
COMT I NlJE 






o n 


// '^U"> 



r> I '-'F vs T O--' L (1 7 ) * A ; 1 7 * ?6 ) » D ( 1 ? ) ♦ F P-n V ( 1 2 ) f F I ^AAX ( 1 2 ) 

= VVAR + ^ 


DC J=1,I 

J<=J + 1 


M=LC^) 

r { D ) R A , 56 * 74 ^ 


D{-vM = 4 ( »^ •'■■’VARD'^FI- IN( y ) 

r,>' TC' ^c. 


QVT ! VUF 
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// OUo 


i^i.f.Tor?r 


S U R p* 0 U T I - soil i'.: ( N V A R 1 » N V A p. » ' M 

LIMFAR PRRRPA^-^a^!^, FOR THE 1130 — SOLUTIO>J PHA5F 

D I MPVC; y rjf.; n(lR) 

CON^^IQNi A n 7 I / A ) t ■•;' ( 1 7 ) ^ L ( 1 7 ) * TI ! » J J ^ II > I W 

covyo'-: D I VAAX ( 1 ? V * T ''M ( 1 ? I » F^C ( 1 P ) *00 ( ] R ) 

A r.^ A X ( 1 2 ) * F T 0 I N ( 1 2 ) * C ( 3 * 1 3 ) * T C ( 3 ) » 0 i/r » T C QfJ^±CQ^) 

7 rOR'-AT ( 7 I R ♦FP . A ) 

3 PpR^AAT ( T 4* F?0*B ) 

A rnpy A"!* ( i i‘ ! i -'FEA s iol^ ) 

7 F';R'V>T < 7 4H \/AoIAPLF V ALUFJ ^ 

« FORMAT ( 1 1 ■ < F'J" C T I ONAL * F 2 0 . B * / / ) 

^ FOR''^A‘’' ( ?AH VA'-' I ARLF SMAO* COST) 


H H L R R 0 G 
V Lr’ K OG 
BOLPROG 
3 1 LPRQG 
'5 2LPROG 
b?LPROG 


=^0R'-'AT(Oh ^'FASIFLE} 
^C^^’AT (IA,^1F.2*10X*IA^ 
F f:R A T { 1 0 LJM 0 ^ n ) 

CALL D A T -S ( 7 . M I SW) 


GO TO ; ♦ 701 ) sv; 

CONTIMUE : 

v.-RI IV, 101 ) 

C 0 N! T I "■! 1 1 ^ . 

'/ Kv'sQ 


0 TO (An,F2)*I5W 


IF(L{ I ) ) 77 *2F,23 
0: 


IF ( A( I * J ) ) 26*27 *2ft 
A n I I * J) =M I I T * J)-AM_ 
COAT Pi! jP 

GO T ^ 7 3 

i^ = I I I 


F6LFR0G 

F 7LPR0G 

FBL^^ROG 


5 F LPRQG 
6 0 L P R 0 G 


62LPROG 


6ALPPC6 


66LPR0G 

67LPR0G 


68LPR0G 
69 LPRQG 
70LPROG 

7 IL PP0 G 

72LPROG 


?4LPR0G 


I F ( A ( X , j } ) A3 *A2 *A2 
T p { ■.. (■ - / ) ^ A , j ] ) A 2 » A 2 * A 7 
a' { ) =A { ■«' * J ) 

^ ly \ — I 


GO TO 4? 

46 IF ( L CO ) 44 *67*A6 

46 F J = L(>^ ) 

^ r— n.Tt 


IF ( A ( t *F J ) ) 170* 120*1?] 

o CO NT IN 


76LPR0G 

77LPRGG 


78LPROG 
79 LPRQG 
80LPROG 
SI LPRQG 


62LPROG 
83LPROG 
84LPR0G 
8 5 LPRQG 


86UPROG 

87LPRWG 


C-19 








Ar,F niA 

i,(p TTFfTM-T'l^ni 




STOP 

' 


89LPPWG 
90LPR0G . 

^ -i- 1 — — 

c. T = T + 1 

’ 


9 ilppog 

92LPHCO 

IF{ I-II )52»'^?*56 
(\ 7 TF/A^T-K'f^>Fri*Sn«*5l 



9 3 L P R 0 G 
^4LPR'.-^G 

^ iC. ir\f\\i%^ui / .' '■■ ■ ■ y •■'V * ■': it... — — 

SI X = A( I *JJ) /At I *KJ) 



9 5 L P R 0 G 
96LPRWG . 

jr _ ;— — 

S5 IFt 53»S0»50 

c, Y f .1 T Tv^ ~ V 



97LPRCG 

98LyR0G 

JK=I 



99LPR0G 

ir-OLPROG 

I 2 J LL !lLL — 

S6 X = A ( J«< .k' J) 
t f 1 \ 1 



lOlLPRCG 

102LPR0G 

L ^ J ^ ^ — 

DO 5 7 I = 1 • I I T 

7 t ^ — A / T fc M 



103LPP0G 
1 04LPRWG 

^ / 'M I } - 1 1 1 — 



105LPROG 
106LPR0G _ 

DO 1 a 1 • 1 vJ 

DO 59 J=1^JJ 
Tc/A/iY-n\Fo»5P*5P 



107LPR0G 

KiFLPROG 

} r? , * ^ y 7 ^ — 

55 IF ( W(I ) ) 5P0f59*5«0 
con AfT. M— A/T« )i *-'«/ (I)-W-tA(JK*%J)/X) 



1 09 L PROG 
irOLPROG 

59 COMTINUE 

T f — 1 Y a. 1 



lilLPROG 

112LPR06 

i J ^ J ^ 1 ' 

DO 61 I=IJ*III 

T V« ) I 



113LPR0G 
114LPR0G _ 

- O - n 1 ^ ^ '. 

I F ( At J‘^ * J ) ) 60, 61 4 60 

A TP/ Ul / TT^AnO*Al«600 



1 15 LPROG 
116LPR0G 

6 0 A ( I , J ) «A ( I * J ) “W ( I )*(A(JK»J)/X) 

A 1 rnMT T Ml !P 



117LPPOG 
liRLPROG ^ 

ox >v. -J ;N t i !N c — — 

DO ?05 J=1 * JJ 

OOP A/ lY * n = A< \^*J)/X 



119LPR0G 

120LPROG 

i< ]<r < Y + l 

nn TO ( 707 • 705 ) I 



121LPR0G 

7'>5 CO^'TINlIH 

MO TTFt TVUI 061!<:KK,A(X,JJ.)*L(J<) 



122LPRCG 

702 COMTIMUE 

r A 1 1 OAT Fiof ( Xi • T SW } 



123LPROG • 

0 L L. / \ 1 D ^ \ T i V »T t 

GO TO (70,^A),ISV' 

4.0' TPfY-.1 T70.7n«63 



124LPRCG 

125LPRCG 

63 IJ=JJ-1 

nr; p p != 1 • T . J 



126LPROG. 

127LPRQQ.' 

IF(A(<,J)-*00nl)65, 65*66 



128LPROG 
1 29LPRQG 

60 rOMT T '1 Ih 

■A'PJlTFt IW*103) 

r\/> 1 — 



130LPR0G 

131LPR0Q 

nC; 1 1 J 1 * Jii - 

15 i A ( I I I * • 0 



132LPROG 

133.L£RQG__. 

^ t - — — 

r r\ T rs /l i± 




66 i/'RlTEt IW*6) 



136LPROG 
137LPROG 

P AU.Ob — — 

70 WRiTEt IW*F )A( 1* JJ) 
r,n TO { 7 04 • 705 ) * 1 5W 



138LPR0G 

7n«; cONTir-UJE 

yg!TE( IW.7) ^ ^ 



139LPROG 
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riQ 71 I = ?4.n 

7 : V."7 ! TE ( I * 5 ) L n ) t A ( I ; J J ) 

4 CQNTIMU^ 

CALL DATS''.' ( I 5W ) 

GQ TO ( 77.7SUISV/ 

17 . W9 I T^~ ( IW ) 

IJ = JJ-1 

no 74 J=1 * IJ 


73 '/R ITE ( IWf S ) J»A( 1 » J ) 

74 CO^:T!MiiF 

75 CALL RATS'^' { ? * I SW ) 

GO TO 

75 TJ=JJ-1 

DUV^V-Y QPAO TO FFFD R1 AN 
REAn(2,5) 


no 77 1=1*111 

IF(A( 1 >Jn 73>77>7_E 

7-" V’R I TE ( R * 7 ) I ♦ J* A ( I * J ) 

11 rovTlNlF 

U'RITC( 2, 100 ) 


RO WR r TF ( 7 * ? ) I »L( I) ( I * JJ) 

0 CALL DAT.S^M At T SW ) ^ 

GO TO (B1 *P2) * ISW 

f?r Ti iPM 

2 CALL OATS'" (5* I SW) 

r-O TO (B3.P4).ISW 

CALL LTN>^(oIGHT) 


14QLPRQG 

141LPR06 


NG STATION THJ 


142LPROG 

I43LPRQG 

144LPR0G 

IASLPR 'OG 

146LPP0G 

1 / -» I r> /"i ^ 


148LPRCG 
l49 LP ROG 
ISOL^ROG 
ISILPR^ G. 
152LPROG 
153LP RGG 
154LPR0G 
■ 1 5 5LP R .Q^ 
1 5 6 L P R W G 


158LPPGG 

,I5S‘LPRC..Ci 

16clppv:g 


162LPR0G 


17CLPRCG 
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Appendix D 


SHUTTLE NAR DATA 



NORTHROP SERVICES* INC. 


TR-1246 


SHUTTLE ASCENT DATA MAX Q 


m = 1,313,800. 

C^-O.. 

= 30,812,000. 


I = 12,682,000. 


I = 316,950,000. 


S ^ = 249.91 
ref 

*')l«e23L “ 

b . = 32.8 
ref 


i = 31,233 

"we4L = -^502 

g = 9.77 

C - „ - -.3165 

£6e3L 

e = 53.2 
0 


U = 398. 
o 

- 0- 

W = 7.55. 
o 

Cng= 1.249 

= 398. 

C' =-.1019 
nor 

Cyg = -2.567 

C . = -.1831 

n6el 

= .1092 

Yor 

See23L = '*377 

Seel " 

S6e4L = 

See23L ^ *^5^2 

C . „ = -1.158 
noe5L 

S6e4L = 

C =0. 
nr 

See5L = 

C =0. 
np 

Sr = 

= .00035 nt/deg 

% 
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NORTHROP SERVICES, INC. 


NAR CONTROL SCHEME MAX Q 


Torque Command 


= ^ip ^ ^ 


i = r + <- + 


otf; 


• Signal Distribution Scheme 
6el = 6e23D = 2<|> 

6R = 5^ 

6e23L =i 


6e4L = V 2 ~ Y 


6e5L = V 2 (ip + 4 (|>) 


D-2 



